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Abstract 

The neutral B meson lifetime has been measured with the data collected by the BABAR detector 
at the PEP-II storage ring during the year 2000 for a total integrated luminosity of 20.3 fb . The 
B° — > D*~ir + decays have been selected with a partial reconstruction method in which only the 
fast pion from the B° decay and the slow pion from D*~ — > D°ir~ are reconstructed. The B° 
lifetime has been measured to be 1.510 ±0.040 ±0.038 ps with a sample of 6971 ±241 reconstructed 
signal events. 
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1 Introduction 



The technique of partial reconstruction of D*~ mesons (charge conjugate states are always implied), 
in which only the slow pion from D*~ — > D tt~ is reconstructed, has been widely used in the past 
H] to select large samples of reconstructed B mesons. This technique provides a way to measure 
the combination of the CKM unitarity triangle angles (2(3 + 7) with B° — > D*~n + decays [Q]. The 
present measurement has been performed as a first step towards the goal of an analysis measuring 
the angle 7. In this respect, the reconstruction of the signal events, the rejection of the background, 
the characterization of the various background components and finally the study of the At resolution 
obtained in these events are important tools for a CP analysis. All these tools are presented here 
in the context of a B° lifetime measurement. 

2 The BABAR detector and dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage 
ring during 2000 and correspond to an integrated luminosity of 20.3 fb _1 collected at the T(4S) 
resonance and 2.6 ftr 1 collected 40 MeV below the resonance for background studies (off-peak 
events) . 

PEP-II is an energy asymmetric storage ring, with positron and electron beam energies of about 
3.11 and 9.0 GeV. The center-of-mass frame of the e + e~ collision is therefore boosted along the 
z direction in the laboratory frame, enabling decay time-dependent measurements of B mesons 
through vertex reconstruction. 

Samples of simulated BB and continuum events were analysed through the same analysis chain 
as the data. The equivalent luminosity of the generic simulated data is approximately equal to one 
third of the on-resonance data, while the equivalent luminosity of a specialized simulated sample 
containing B° — > D*~ir + decays followed by D*~ — > D°ir~ is 2.9 times larger than the on-resonance 
data sample. 

A detailed description of the BABAR detector and the algorithms used for the track reconstruc- 
tion, particle identification and selection of BB events is provided elsewhere Q; a brief summary 
is given here. 

Only charged particles are used for the partial reconstruction of the signal. Particles with mo- 
mentum higher than 170 MeV/c are reconstructed by matching hits in the silicon vertex tracker 
(SVT) with track elements in the drift chamber (DCH). Since tracks with momentum below 
170 MeV/c do not leave signals on many wires in the DCH due to the bending induced by the 
magnetic field, they are reconstructed in the SVT alone. 

Electron and muon identification is used in veto mode for the selection of the fast pion. Electrons 
are identified on the basis of the energy deposited in the electromagnetic calorimeter (EMC), the 
track momentum and the energy loss in the DCH. Muons are selected by requiring deep penetration 
in the instrumented flux return (IFR). 

The Cherenkov light emission measured in the particle identification detector (DIRC) is em- 
ployed to reject kaons from the fast pion sample. 

Neutral particles are reconstructed from clusters in the EMC that are unmatched to projected 
charged tracks. They are used only to compute event shape quantities. 
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3 The partial reconstruction technique 



The B° — ► D*~ir + decays are reconstructed using only the tracks from the 7r + (fast pion) and from 
the 7T~ (slow pion) in D*~ — > D°7r~ . The fast pion momentum in the Y(4S) rest frame is required 
to be in the range kinematically allowed for the decay B° — > D*~ir + . For nominal values of the 
beam energies this momentum is in the range 2.114 - 2.404 GeV/c. The momentum of the slow 
pion is required to be greater than 50 MeV/c. 

Assuming that the fast and slow pion come from the decay B° — > D*~ir + followed by the two 
body decay D*~ — > Xir~ , it is possible to compute the mass of the recoiling system X averaged 
over an unmeasured angle due to the unknown direction of the B momentum in the T(4S) rest 
frame. For B° — > D*~ir + decays this recoil mass M rec peaks at the D° mass with a width slightly 
less than 3 MeV/c 2 . Assuming that the B momentum lies in the plane defined by the fast and slow 
pions momenta in the Y(AS) rest frame, it is possible to compute the helicity angle of the pion in 
the D*~ rest frame as well as the candidate D direction. 

Since the dominant source of background is continuum events, the selection procedure aims at 
reducing this contribution. The main requirements of this selection are the following: 

• R2, the ratio of the second to the zeroth Fox- Wolfram moment [Q], computed from charged 
particles, is required to be less than 0.35. 

• No other tracks should be in a cone of opening angle 0.4 rad centered on the fast pion 
momentum in the T(4S) rest frame. This cut is effective against continuum events because 
in this case tracks tend to be clustered in jets. 

• A Fisher discriminant Fd is computed from 15 event shape variables. Among these are 
the scalar sum of the momenta in the T(4S) rest frame of all tracks and neutrals, in nine 
20° angular bins around the fast pion direction. The signal peaks at Fd = —0.5 while the 
continuum background peaks at Fd = 0. Events are required to satisfy Fd < —0.1. 

• The cosine of the soft pion helicity angle is required to be larger than 0.4 in absolute value. 

The B° decay point is determined from a vertex fit of the fast and slow pion tracks and the 
beam spot position in the plane perpendicular to the beam axis (the x-y plane). The beam spot is 
determined on a run- by-run basis using two-prong events Q. Its size in the horizontal direction is 
120 /zm. Although the beam spot size in the y direction is only a few ^m, a beam spot constraint 
of 30 /jm is applied to account for the flight of the B° in the y direction. Only events for which the 
probability of the vertex fit is greater than 0.1% are considered further. 

The decay point of the other B is determined from a subset of the remaining tracks in the 
event. All the tracks with a center-of-mass angle greater than 1 rad with respect to the candidate 
D° direction are considered. This requirement is used to remove most of the tracks from the decay 
of the D° daughter of the D*~ , which would otherwise bias the reconstruction of the other /B/ 
vertex position. The selected tracks are then constrained to the beam spot in the x-y plane. The 
track with the largest contribution to the vertex % 2 , if greater than 6, is removed and the fit iterated 
until no tracks fails this requirement. Vertices composed of just one track and the beam spot are 
rejected in order to reduce the number of poorly measured vertices. Simulation shows that after 
all these requirements in about 85% of signal events the other vertex has no tracks from the D° 
decay. 

Figure |l] shows the recoil mass distribution for on-resonance data events obtained by applying 
all the above cuts when the fast and slow pions have opposite-sign charges. These events are used to 
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measure the B° lifetime. Events with same-sign charges for the two pions are used as a background 
control sample (same-sign sample). Another background control sample is obtained by reversing 
the cut on the Fisher discriminant and requiring > (BB depleted sample). 

The individual distributions shown in Figure [j] are obtained by fitting simultaneously to the 
on- and off-resonance data the contributions from 

• B° -> D*-ir+ events; 

• B° -> D*-p + events; 

• B° -» D*-£ + u events; 

• peaking BB background, where the fast and slow pions originate from the same B meson 
decay; 

• non-peaking BB background, excluding the events from the previous four categories; 

• peaking cc background, where the slow pion comes from the D*~ — » D°ir~ decay; 

• continuum background excluding the events from the previous category. 
All peaking backgrounds peak in the recoil mass distribution. 

The shape of each contribution has been determined from the Monte Carlo simulation while all 
the normalizations have been left free in the fit except those for B° — > D*~p + and B° — > D*~l + v 
events which have been fixed to the value from the Monte Carlo simulation. The B — ► D*~ p + 
branching fraction is taken to be (6.8 db 3.4) 10 -3 || giving an estimate of 1570 events in the final 
sample from this decay mode. The fit gives 6971 ±241 B° — > D*~ir + events. In the lifetime analysis 
to be described later, the B° — ► D*~ p + events will be included as signal in addition to the 6971 
B° — ► D*~ir + events. Table [l] reports the composition of the selected sample for recoil mass larger 
than 1.86 GeV/c 2 (signal region sample) as determined by this fit. 

Table 1: Composition of the data sample in the signal region. 



Source 


Fraction in Signal region (%) 


£ u -► D*~TT+ 


46.1 


B° -» D*~p+ 


9.3 


B u -» D*~l + v 


2.3 


BB peaking 


8.9 


BB non peaking 


11.5 


uds 


12.7 


cc non peak. 


5.2 


cc peaking 


4.0 



4 The lifetime measurement 

The PEP-II collider produces BB pairs moving along the beam direction {z axis) with an average 
Lorentz boost of (/?7) = 0.55. The lifetime is determined by measuring the quantity — Zdecay — 
Mother, where Zdecay (z ther) is the position along the beam line of the reconstructed B — ► D*~tt + 
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Figure 1: Recoil mass distribution obtained with the selection explained in the text when the 
fast and slow pion have opposite-sign charges. In this plot the small B° — ► D*~l + u component 
has been added to the peaking BB contribution. The continuum background includes the peaking 
cc component. The fit gives 6971 ± 241 B° — > D*~tt + events. The region to the right of the line 
corresponds to the signal region. 

decay (other) vertex. To remove badly reconstructed vertices, all events for which a^ z > 400/Ltm, 
where a& z is the uncertainty on Az computed for each event, are rejected. 

Residual tracks from the D° decay, not removed by the track selection for the other B vertex, 
bias the reconstruction of the other B vertex position. This bias is removed by correcting the Az 
value for each event with a correction function determined from the simulated signal sample as a 
function of Az. 

The proper time difference between B decays is then computed with the relation At = Az/{cf3~f). 
A fit with a double Gaussian to the At residuals in the Monte Carlo simulation shows that 75% of 
the events are contained in the narrower Gaussian, which has a width of 0.8 ps. 

The lifetime t b o is obtained from an unbinned maximum likelihood fit to the two-dimensional 
At, a At distribution. The At distribution of signal events is described by the convolution of the 
decay probability distribution 

f(At true \T B o) = _L e -|AW|/r B o j (!) 
2t b o 

with the experimental resolution function, which is represented by the sum of three Gaussian 
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distributions. The first two, accounting for more than 99% of the events, have the form 



G(S(At), a At ) = 1 e ^T" (2) 



where 5 (At) = At — Att rU e is the difference between the measured and the true value of At, b is a 
bias due to the charm tracks in the other vertex and the scale factor S is introduced to account for 
possible misestimation of the error on the proper time difference At. The third Gaussian of fixed 
bias 6 = and scale factor S = 6 is added to account for badly mismeasured events ( "outliers" ) . 

The B° — > D*~p + resolution function has been found with Monte Carlo simulation to be the 
same as for signal events. In the fit for the B° lifetime the B° — > D*~ p + events are considered 
as signal events. Events coming from the decay B° — ► D*~l + v are added to the BB peaking 
component and the At distribution of the latter is assumed to be equal to the BB non-peaking 
component. 

The At distribution of the B background contribution is described by a resolution function 
similar to that used for the signal. In this case the lifetime parameter represents an effective B 
lifetime which has been fitted independently from the signal lifetime. 

The At distribution for uds events is described by the sum of three Gaussians while the At 
distribution for cc events has been described by the sum of a Gaussian plus the same Gaussian 
convoluted with an exponential term to account for the effective charm hadron lifetime. 

The function used to fit the data is the weighted sum of four contributions: 

F(At,a,T B o) = [1 - f BB (M rec ) - fuds(M rec ) - f c5 (M rec )]F B o(At,a,T B o) + 
f BB (M rec )F BB (At,a) + f uds (M rec )F uds (At,a) + f c5 F c5 (At,a) 

where the functions F B o, F BB , F u d s and F c£ describe the measured decay time difference distribu- 
tions for the signal, BB background, uds and cc events, respectively. f BB , f u ds an d fee are the 
probabilities that the event is from the BB, uds and cc background, computed for each event on 
the basis of the measured value of the recoil mass M rec . 

The key parameters describing the At distributions for the background events are fitted on the 
BB depleted control sample for the continuum background and on the same-sign sample for the 
BB background. The BB depleted control sample has been used because the integrated luminosity 
of the off resonance sample in not sufficient to precisely determine the parameters of the continuum 
background. An alternative parameterization obtained by fitting the side-band control sample has 
also been tried giving compatible results. 

All the parameters describing the signal resolution function are free in the fit except the bias 
of the second Gaussian which is fixed to the value found on the signal Monte Carlo sample. 

The result of the fit to the signal region sample in the range |At| < 15 ps is t^q w = 1.524 ± 
0.040 ps, where the error is statistical only. Figure 2 shows the comparison between the measured 
At distribution and the fit result. All the fitted parameters appear in Table 

The raw value r B % w value is corrected for a small bias, 0.014 ± 0.020 ps, observed when fitting 
the signal Monte Carlo sample using the same procedure, yielding the corrected value t b o = 
1.510 ±0.040 ps. 



5 Systematic uncertainties and cross-checks 

The systematic error on t b o is computed by adding in quadrature the contribution from several 
sources, described below and summarised in Table ||. 
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Figure 2: Decay time difference (At) distribution for the signal region sample on a linear (upper) 
and logarithmic (lower) scale. The curve shows the result of the unbinned maximum likelihood 
fit. The shaded areas represent from bottom to top the following contributions in the fit: uds+cc, 
non-peaking BB, peaking BB and signal. 



The fractions of the various background components are varied by their uncertainties obtained 
from the recoil mass fit. The parameters of the various background At distributions are also varied 
by their uncertainties, properly accounting for their mutual correlations. The effective lifetime of 
the BB peaking component is varied by ±0.044 ps corresponding to the difference between the 
value found in the data and in the Monte Carlo sample. An alternate resolution function for this 
component has also been tried with negligible variation on t b q. 

The bias of the wide Gaussian, the only parameter of the signal resolution function which is 
not fitted, is varied in a conservative range. Several different analytical expressions are used to 
represent the small fraction of outliers. The fit range is varied from (—10, 10) ps to (—20, 20) ps. 
The parameters of the Az correction are varied according to the uncertainty due to the finite signal 
Monte Carlo sample size. The systematic uncertainty related to this correction is estimated from 
the effect of a ±5% variation on the fraction of fitted other- i? vertices which have no tracks from 
the D° decay. 

The z length scale is determined to about 0.4% from secondary interactions with a beam pipe 
section of known length. The statistical uncertainty of the residual bias found on signal Monte 
Carlo events is also added to the systematic error. 

The total systematic error of ±0.038 ps is found by adding in quadrature the uncertainties from 
the above sources. 
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Table 2: Result of the fit to the data. The uds and cc backgrounds have been fitted on the BB 
depleted control sample. The BB background has been fitted on the same-charge control sample. 



Parameter description 


Value 


uds narrow Gaussian fraction 
uds narrow Gaussian bias 
uds narrow Gaussian scale 
uds outliers fraction 
cc narrow Gaussian bias 
cc narrow Gaussian scale 
cc lifetime fraction 
cc lifetime 


0.972 ± 0.079 
0.181 ±0.020 
1.194 ±0.037 

0.0070 ± 0.0028 
0.342 ± 0.039 
1.404 ±0.091 
0.296 ±0.102 

0.784 ±0.196 ps 


BB lifetime 

BB narrow Gaussian fraction 
BB narrow Gaussian bias 
BB narrow Gaussian scale 
BB outliers fraction 


1.611 ±0.044 ps 
0.721 ±0.104 
0.326 ±0.131 
0.863 ±0.210 
0.000 ± 0.013 


raw £?° lifetime 

signal narrow Gaussian fraction 
signal narrow Gaussian bias 
signal narrow Gaussian scale 
signal wide Gaussian bias 
signal wide Gaussian scale 
signal outliers fraction 


1.524 ±0.040 ps 
0.962 ± 0.065 
0.140 ± 0.058 
1.336 ± 0.082 
-0.5 (Fixed) 
2.640 ± 1.420 

0.0049 ± 0.0078 



The dependence of the result on several different variables (angular width of the cone used 
to reject the D° tracks, fast pion momentum, polar and azimuthal angles) has been studied: no 
statistically significant effect is found. 

6 Conclusion 

In conclusion the neutral B meson lifetime has been measured with a sample of 6971 ± 241 partially 
reconstructed B° — ► D*~ir + decays: 

t = 1.510 ± 0.040 ± 0.038 ps. 

This preliminary value is consistent with other recent BABAR measurements || and with the world 
average B° lifetime [||. 
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Table 3: Summary of contributions to the systematic error. 



Source 


Error (fs) 


Background Parameters 


13 


Fractional composition 


20 


Az correction (MC) 


6 


Az correction (model) 


17 


Bias of the wide Gaussian (signal) 


6 


Outliers 


9 


At range 


9 


MC bias 


20 


z scale 


6 


Total 


38 
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